arXiv:1501.00806vl [nucl-th] 5 Jan 2015 


Analytical solution for the Davydov-Chaban Hamiltonian 
with sextic potential for 7 = 30° 

P. Buganu^’Q and R. Budaca^ 

^Department of Theoretical Physics, 

National Institute for Physics and Nuclear Engineering, 

Str. Reactorului 30, RO-077125, POB-MG6, Bucharest-Mdgurele, Romania 

(Dated: January 6, 2015) 

Abstract 

An analytical solution for the Davydov-Chaban Hamiltonian with a sextic oscillator potential 
for the variable /3 and 7 fixed to 30°, is proposed. The model is conventionally called Z(4)-Sextic. 
For the considered potential shapes the solution is exact for the ground and {3 bands, while for the 
7 band an approximation is adopted. Due to the scaling property of the problem the energy and 
B{E2) transition ratios depend on a single parameter apart from an integer number which limits 
the number of allowed states. For certain constraints imposed on the free parameter, which lead to 
simpler special potentials, the energy and B{E2) transition ratios are parameter independent. The 
energy spectra of the ground and first ft and 7 bands as well as the corresponding transitions, 

determined with Z(4)-Sextic, are studied as function of the free parameter and presented in detail 
for the special cases. Numerical applications are done for the 128,130,132^^ I92,i94,i96p^ isotopes, 

revealing a qualitative agreement with experiment and a phase transition in Xe isotopes. 
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I. INTRODUCTION 


Soon after the Bohr-Mottelson Model (BMM) [l|, l2|] was proposed for nuclear structure 


together with its hrst solution |l| for spherical nuclei, many attempts were done to improve 
and extend it by taking into account axial and non axial deformation, coupling between 
/5 and 7 vibrations or various anharmonicities. Most of these approaches were reviewed 
in Refs.js-Sj. A new phase in the field begun with the studies of the phase transitions 


III21 


by means of the classi^! 
compact Lie algebras 


limits of Hamiltonians constructed from operators belonging to 
8 |. This algorithm was further used in {q] to associate classical 


shape variables to the Interacting Boson Model (IBM) 10l-ll3l| and allowed to establish the 


nature of the quantum phase transition 


14| between the dynamical symmetries, namely. 


U(5) (spherical vibrator), 0(6) ( 7 -unstable) and SU(3) (axial rotor). The start of a long 
series of studies, both theoretical and experimental, was given mainly by two papers in which 
approximate solutions of the BMM were offered for the critical points of the shape phase 
transitions U(5)-0(6) and U(5)-SU(3), called E(5) jl5| and Xf 5l llT l. respectively. Other 
two important models which are worth mentioning here are Y(5) [l7| and Z(5) 18| associated 
with the transitions between the axial and triaxial shapes and respectively between prolate 
and oblate shapes. The critical point approaches mentioned above have the advantage 
to be parameter free solutions except for a scaling factor, making them easily verifiable 
reference points for the experimental data. This is actually a general trait of the exactly 


solvable models of nuclei 


19l |. Other efforts were also directed to special realisations of the 


BMM in the view of some constraints imposed on the shape variables or inertial parameters. 
For example ’’freezing” the 7 variable to a certain value in the classical BMM, leads after 
quantization in curvilinear coordinates to simpler Hamiltonians suitable to describe the 
special case of the 7 rigid collective motion. The hrst study in this direction brought to 
light the Davydov-Chaban model for rotation-vibration interaction in non-axial nuclei 20|. 
Later on, an exact solution for this model was proposed 2l|] in the case of 7 = 30°, where 
instead of a displaced harmonic oscillator in (3 shape variable an inhnite square well potential 
was used. The solution called Z(4) due to the similarity to the Z(5) model, inspired other 
studies of the 7 rigid solutions 2^24|. 

In this paper we propose an analytical solution for the Davydov-Chaban Hamiltonian 
with 7 = 30° and a sextic potential for the only shape variable, i.e. f3. The model is 
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conventionally called Z(4)-Sextic. In this framework, the separation of the angular variables 
from the /3 shape variable is exact. The differential equation involving Euler angles is solved 
as in Ref. 2^, while that for /3 is brought to a Schrodinger form with a sextic potential 


and a centrifugal-like term. The problem of the sextic potential is not an exactly solvable 
one because its spectral problem is reduced to the diagonalization of an inhnite-dimensional 
Hamiltonian matrix. However, for a family of potentials whose coefficients satisfy certain 
relations between them and the factor of the centrifugal term, the problem becomes quasi- 
exactly solvable j3, 27|, i.e. its inhnite Hamiltonian matrix acquire a block diagonal form 
allowing thus an algebraic treatment for a hnite subset of eigenstates. For a physically 
meaning description, the above mentioned constraints must be corroborated also with the 
condition of constant potential. Despite these restrictions, the Z(4)-Sextic eigenvalue prob¬ 
lem is exactly solved for the ground and (3 bands. Concerning the 7 band, an approximation 
is involved in the centrifugal term in order to accommodate all the above restrictions. Due 
to the scaling property of the exactly solvable sextic potential with an associated centrifugal 
term, the energy and the B{E2) transitions depend on a single parameter up to an overall 
scaling factor. Moreover, for particularly interesting shapes of the potential, parameter free 
expressions are possible for the normalized energies and B{E2) transition probabilities. 


The use of such an involved potential is supported by the fact that it is the simplest shape 
which through continuous variation of its parameters can have either a spherical minimum, 
a deformed minimum or both. It is worth to mention that exact 28|, [2^ and approximate 
3, M| solutions by using a sextic potential were also given, in Eve dimensions, for E(5) and 


respectively X(5) and Z(5) related approac 


les. Other solutions in the vicinity of 7 = 30°, 


but with 7 soft can be found in Refs. 32h36|. 


The present work has the following plan. The Z(4)-Sextic model Hamiltonian is presented 
in Section H, while its associated (3 differential equation is treated in Section HI. In Section 
IV, one gives the model wave functions and calculate the B{E2) transition probabilities. 
Extensive numerical results and few model hts to experimental data are given in Section V. 
The main conclusions are drawn in Section VI. 
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II. THE MODEL HAMILTONIAN 


The eigenvalue problem, when the nucleus is 7 —rigid, has the following form 20| : 


_L V 

2B [/33 3(5 4/32 


vI>(/?,0) + f^(/3)vI>(/?,0) = ET(/3,0), (2.1) 


where B is the mass parameter, /3, 7 and Qk are the intrinsic deformation coordinates 
and respectively the operators of the total angular momentum projections in the intrinsic 
reference frame, while with O are denoted the rotation Euler angles ( 01 , 02 ,^ 3 )- Here, 7 
is considered a parameter and not a variable, such that when the kinetic energy of the 
classical BMM is quantized in curvilinear coordinates one arrives at the Hamiltonian (12.Ih 
which depends only on four variables (/5, H). When 7 = vr/b, two moments of inertia in the 
intrinsic reference frame become equal and then the rotational term reads: 

Q\ .Ao 3 


1 3 
4 sin^ 


= (Q -yQf)- 


( 2 . 2 ) 


The separation of variables is achieved by considering the wave function 4/(/5, hi) = (j){/3)'ip{Q) 
which leads to the following equation in /3 variable: 


I d ^3 d W 
' 13^ dl3^ dfd ^ /32 


0(/?) = e0(/3). 


(2.3) 


where the following notations are used n(/5) = ^V{(3) and e = while W is the 

eigenvalue for the equation of the angular part, 
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The above equation was solved in Ref.[25j with the results: 
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W = Wlr = L{L + 1 ) 




and 


V’i'K(H) = 


2L + 1 




( 1 ) 


(2.4) 


(2.6) 


( 2 . 6 ) 


167r2(l + 6 r ^ o ) L 

Here £*7(^) the Wigner funetions associated to the total angular momentum L and 

its projections on the body hxed x-axis and laboratory hxed z-axis, R and respectively /i. 
For the energy spectrum it is more advantageous to use instead of R the wobbling quantum 
number = L — R which for the ground and fd bands is = 0, while for the 7 band 
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it takes the values = 1 for L odd and = 2 for L even. Within this convention the 
eigenvalue of the angular part of the problem is written as 


»LH = Wl„. = L(L + 1) - -(L - n„f. 


(2.7) 


III. SOLUTION FOR THE j3 PART OF THE HAMILTONIAN 


It is convenient to bring Eq. fl2.3D to a Schrodinger form. This is realized by changing 
the function with (j){/3) = j3~iip{(3) 120 1: 




+ 


+ v{^) 




Further, Eq. fld.Sp is compared with the exactly solvable case of the sextic potential 
which leads to the following correspondences: 


(3.8) 


Wlu^ 


= 2s 


\)i-l 


v{f3) = 


Aa i s + - + M 


2 06 


( 3 ^ + 2a6/3^ + a"/? 


(3.9) 


(3.10) 


The potential fl3.10p depends on two parameters, a and 6 , and on L and quantum 
numbers through s. M is a natural number which establishes the number of states that can 
be determined. This implication will be explained later when discussing the wave functions. 
The number of parameters is reduced to a single one by changing the variable with {3 = ya~^. 
Then by introducing the notations a = -^ and £y = -^ one gets: 


d^ . Wr^ + 


dr 


+ 


y" 


+ — Ac)y‘^ + 2ay^ + y^ 


v{y) = £yV{y), 


where 


c = s 


1 

2 


M. 


(3.11) 


(3.12) 


Because s depends on L and the potential of Eq. fl3.1ip is state dependent. For the 
ground and (3 bands = 0, such that 


^3 L 3 

*= 4 + 1 ’ * = "+ 4 + 2 - 


(3.13) 


In order to have a state invariant potential for this case of ground and /3 band states, the 
following condition must be satished: 


L 3 

c = M + — + - = const. 
4 2 


(3.14) 
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It can be easily checked that the above condition is satished if M is decreased with one nnit 
when L is increased with fonr. This means that for L/2 even and L/2 odd there are two 
different constants c: 

(M, L) : {K, 0), (i^ - 1,4),^ + ^ = c^, (3.15) 

(M, L) : (i^, 2), (ir - 1, 6 ),... ^ + 2 = cf, (3.16) 

which differ from each other jnst by 1/2. Note that the valne of = Mmax pnts a limit 

on the nnmber of states which might be determined. For example if iF = 1, the maximnm 

angnlar momentnm state which conld be analytically described wonld be the L = Q state, 
while for K = 2, the L = 10 state and so on. This is actnally a direct conseqnence of the 
condition (I3.14p . In case of the 7 band, when riy^ = 1 and 2, s becomes irrational snch that 
the Eq. (I3.1ip cannot be solved anymore for M integer and with constant potential condition 
fnlhlled. A possible way to handle this problem is to extract from the centrifngal term the 
qnantities 3(L — l/2)/2y‘^ and 3{L — l)/y‘^ for L odd and respectively L even, and to replace 
y"^ with its average {y"^) on r]{y) eigenstates of the remaining Hamiltonian for each angular 
momentum L. With these approximations, s and c from the 7 band will have the same 
expressions fl3.13p as in the case of ground and (3 bands. Moreover, for L even states of the 
7 band, c will have the same two values as for the ground and (3 bands. While for L odd, c 
will have other two values: 

(M, L):{K,l),{K-l,5),...^K + ^ = cf, (3.17) 

(M, L):{K,3),{K-l,7),...^K + ^ = cf, (3.18) 

corresponding to (L — l)/2 even and respectively (L — l)/2 odd. Finally, the four values of 
the constant c can be summarized by the formula: 

cf = cf + = A:+^ + ^m, m = 0,1,2,3. (3.19) 

The condition of the constant potential is then exactly satished for four distinct sets of states, 
which correspond to slightly different potentials. This picture is improved by considering 
for the general potential the following form: 

v^iy) = - 4cf )|/^ + 2ay^ + y^ + wf («)• (3.20) 
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For a fixed K, are constants depending on a, which are fixed snch that the 


minimnm 


energy of the potentials to be the same. Choosing Uq = 0, the other constants are given 
by: 

- 4c*) - (a^ - 4cf) + 2a (i/Jo)" - ( 1 /^:)“ + (j/fo)' - (j/Jj)', 

(3.21) 

with z = 1, 2, 3 and being the minimnm points: 


(</£„)" = k-2a ± ^a^ + 12c*). 


Taking the ansatz fnnction 27(] 


VMiy) ~ PnivV "e "4 1 , 


(3.22) 


(3.23) 


the Eq. fl3.11l) is then rednced to the following differential eqnation, 
4s — 1 d 




d old 
+ 2ay— + 2y^ y— 
dy \ dy 


2M 


PM(r) = APM(r), (3.24) 


y dy, 

where Puiy"^) are polynomials in of order M. The eigenvalnes A are obtained for each 
M nsing the analytical procednre given in Appendix of Ref. Q . For each valne of M there 
are M + 1 solntions which are differentiated by the (3 vibrational qnantnm nnmber np in 
the following way: The lowest eigenvalne A corresponds io np = 0, while the highest to 
np = M + 1. For the present physical problem only the solntions with = 0 and = 1 
will be considered, which correspond to the gronnd and 7 bands and respectively to the (3 
band states. A also depends on L throngh s and one mnst remind that at this point L and 
M are interdependent throngh the condition fl3.14p . the actnal relationship being dictated 
by the valne of K. Thns, the M indexing of A will be replaced from here by K. Following 
all the algebraic manipnlations which lead to Eq. (13.241) and taking into acconnt the above 
considerations, A can be alternatively expressed as: 


A = = £y-u^- 4as - 




{y^)KL {y^)KL 

From the above relation one hnally extracts the total energy of the system: 


3(L-1) 


^nui.2- 


(3.25) 




h\/d 


TLijj ,Z/ 




2B 


3 

A^^l + 4a ( — + 1 ] + + 


F-j). ,3(L-1) 

Onu,,l + 


(y^) 


KL 


(y^) 


KL 


nu,,2 


(3.26) 
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which is indexed by the f5 vibration and wobbling quantum numbers and as well 
as by the intrinsic angular momentum L. The index m is completely determined only by 
L. Although the above energy also depends on the integer K, this is not a true quantum 
number but rather a special kind of parameter. Similarly to the eigenvalue A, the associated 
eigenfunctions of Eq. 03.241) also depend on K and L. Such that due to the orthogonality 
of the angular wave functions 02 .6p . the average of ip' entering in the definition of the total 
energy are only K and L dependent. From Eq. 03.26|) one can see that the energy spectrum 
normalized to the energy of the first excited state depends only on the parameter a and the 
integer K. For further calculations one defines the energy ratios: 


for a fixed value of K. 




Eni3,nu,,L — -^ 0 , 0,0 

-^0,0,2 — -Eo,o,0 


(3.27) 


IV. TOTAL WAVE FUNCTIONS AND B{E2) TRANSITION RATES 


As was explained in Section II, the total wave function is factorized into an angular part 
and a (3 depending factor function: 






(4.28) 


where the angular factor state was defined by Eq. fl2.6|) keeping the notation with R instead 
of Uw for convenience in calculating angular matrix elements. In what concerns the (3 wave 
function, it has the following form: 

= Va(j)npLiy) = VaNj^^L{a)PMnp{y'^)y^e-'T-^ ^ (4.29) 

with y = I3a}/^ and being the normalization constant with respect to the y^dy 

integration measure. As was already mentioned when the expression of the total energy was 
discussed, M is uniquely determined by L for a fixed value of K. Thus, a more natural 
dependence of the total wave function would be on K instead of M. However M express 
more intuitively the analytical form of the (3 factor state. 

Having the analytical expression of the total wave function, one can readily compute the 


B{E2) transition probabili 


as for the Z(4) solution [2l|, 


ies. The quadrupole operator for Z(4)-Sextic has the same form 




(4.30) 
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FIG. 1. The shapes of the energy potential for 7 fixed at vr/h and with K = 1, corresponding to 
a = —10, ±2 y/co(±3.22), 0, 10 are plotted as function of y. 


The reduced E2 transition probabilities are defined as: 




(^2)11,T.M 




(4.31) 


where the Rose’s convention 


37l | was used for the reduced matrix elements. The matrix 


elements over the /3 can be rewritten in terms of y with the following result: 




iM 


M 


iM 


(4.32) 


V. NUMERICAL RESULTS 

As its construction is suggesting, the Z(4)-Sextic model, introduced in the previous sec¬ 
tions, is adequate for the description of triaxial nuclei having a 7 rigidity of 30°. The model 
depends on a single parameter a, apart from a scaling factor and the integer number K which 
gives the extension of the exactly solvable subspace. Depending on the free parameter a 
and regardless of the K value, the sextic potential fl3.20p may have a spherical minimum 
{a > a deformed one {—2^J^ < a < 2^J^) and simultaneously spherical and 

deformed minima (a < —2Jc^) for all considered states. These situations are depicted in 
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Fig. [T] for iF = 1 where one also showed the particular cases of a = 0, a = and 

a = with the latter one corresponding to a potential shape close to that of the inhnite 
square well. The different constants used for the four distinct sets of states depending on 
the parity of L /2 or (L — l)/ 2 , define some small extension intervals where few of the above 
solutions coexists for different states. If all the constants would be equal, the coexisting 
intervals would shrink to a single point value of a. 

The advantage of the present model’s dependence on a single parameter is that one can 
study how its characteristics are changed between the pictures discussed above by contin¬ 
uously varying the free parameter. In order to do this and cover all the above mentioned 
cases, the energy ratios fl3.27p and few interband and intraband B{E2) transitions fl4.3ip 
normalized to the transition 2^ —)■ 0^ are presented in Fig. [2] and Fig. [3l respectively, for a 
sufficiently large interval of a in order to achieve convergence at both sides. The numerical 
results visualized in Figs. [2] and [3] are performed for iF = 1, 2, 3 and 4. For each K there is a 
limited number of available states which are exactly determined in the present model. The 
number of such states in the ground, /3, 7 with L even and 7 with L odd bands increases 
with two when K is increased with one unit. The common parts of the energy spectra cor¬ 
responding to different K are very similar. However, there are some clear differences, such 
as the width and the position of the coexistence intervals identihed in Fig. [2] by the gray 
area where for a set of states the potential shape has a spherical minimum while for another 
set it has a deformed minimum. For the ground and (3 bands the interval is [2 2, 

while for the 7 band the interval is bigger Indeed, as K increases, the grey 

band becomes thinner and its position moves to higher values of a. As a matter of fact 
in this existence interval, one observes a kink in the energy curves which happens at a 
critical value ac- This value corresponds to the absolute maximum of the signature ratio 
-R4/2(«) = i?(0,0,4,Q!) and is interpreted as the critical point for a hrst order shape phase 
transition between spherical and deformed shapes in the framework of presently adopted 
sextic potential. This is also supported by the fact that the hrst derivative of the energy in 
the critical point ac has a discontinuity. As was explained in Ref. j^, the critical point for 
a hrst order phase transition corresponds to the situation when the spherical and deformed 
minima of the potential energy are degenerated. This happens in our case at ac where the 
potential shape is the hattest, being extended over a wide range of non zero deformations. 


A similar critical point was pointed out in the analysis made in Ref. 


29| regarding a 7 
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FIG. 2. The energy spectrum given by Eq. (I3.27h is shown as function of a in the interval [—10,10] 
for iF = 1, 2, 3 and 4. In the panels (a), (c), (e) and (g) are plotted the energy curves of the ground 
band and /3 band which go to infinity when a —oo. While in (b), (d), (f) and (h) panels are 
those corresponding to the 7 band, with the continuous and dashed curves representing L even 
and L odd states, respectively. 
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FIG. 3. The B{E2) transitions 4+ —>■ 2+, 2^ ^ 0^, 0^ 2+ (a) and 3;| —>■ 2+, 2+ —>■ 2+, 
3+ ^ 4+, 4+ ^ 2+, 2^ ^ 2+ (b), normalized to B{E2,2'^ plotted as function of a in 

the same interval [—10,10]. The continuous, dashed, dot-dashed and dotted lines correspond to 
K = 1,2,3 and 4, respectively. 


unstable realization of the sextic potential and in Ref. 38| where a quartic potential was 


involved. The fact that the critical value ac separates two distinct shape phases can also 
be seen from the behaviour of the wave function. For example from the dependence of the 
normalized wave function on a and y, shown in Fig. IH one can see that up to the 

critical value ac the position of the peak changes very rapidly, while for a > ac its position 
do not vary signihcantly. This suggests that the two shape phases delimited by ac have 
different properties. 

Contrary to the energy spectra, the B{E2) transition probabilities shown in Fig. |3]have 
a smooth behaviour as function of a. While the K variation, induce only a small shift to the 
right of the curves from Fig. |3l The common feature of the all considered transitions is that 
their corresponding probabilities become " K degenerate” for a —)■ ±cxd and more sooner for 
the interband transitions. 

As was mentioned before, in the coexistence region, and especially at the critical value 
ttc, the shape of the potential approximated by n ~ 2acy'^ + y^ is the flattest one, which is 
consistent with critical point behavior. Moreover, the potential at ac simulates quite well 
an inhnite square well, supported also by the fact that the corresponding energy spectrum 
is very close to that of Z(4) model. Another interesting aspect of the present model is that 
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0.0 0.5 1.0 1.5 2.0 

y 

FIG. 4. The plot of <j)QQ{y)y^/'^ as function of y for a E [0,10]. The square of this quantity is the 
probability density for the ground state deformation with respect to the dy integration measure. 
The horizontal red line marks the the wave function behaviour at the critical value ac- The 
evolution of the peak’s position with a is visualized by another red curve. The difference between 
two consecutive contour lines is 0.1. 

some energy ratios curves of the ground and fd bands are intersecting each other for a = 0 
and becoming thus degenerate. This can be seen only starting from K = 2, where the last 
two ground band states LMax and LMax — 2 are degenerate with LMax — 8 and LMax — 10 
from the [3 band when a = 0. This degeneracy may reveal some symmetry properties 
associated with the resulting simple potential shape ~ + y^. The low lying 

energy spectrum with a complete set of E2 transition probabilities for this special case is 
graphically represented in Fig. Ofor each considered value of K. A similar representation 
is provided in Fig. [H] for the other special case corresponding to Oc where one also given 
its numerical value. The parameter free results presented in Figs. O and E] can be used in 
a hrst step as reference points for hnding candidate nuclei and then to vary a for a better 
agreement with the experimental data. 
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Another import am 
of the Z(4)-/?^ model 


touchstone of the present formalism represents the exact reproduction 


39l | spectra when a —)■ cxd. This means that the ground and /3 bands are 


degenerated and have a harmonic oscillator type spectrum, while the even and odd angular 
momentum states of the 7 band deviate from this behavior. In what concerns the other 
limit, a ^ — 00 , the ground and 7 band spectra achieve a convergence at a noncollective 
value of i? 4/2 < 2. While the {3 band energy curves go to infinity. The limiting value i? 4/2 = 2 
is reached at a = —0.964, —0.804, —0.666 and —0.545 for A' = 1, 2, 3 and 4 respectively. 

A particular signature which is often used as a characteristic of structure and its evolution 


is the staggering in the 7 band energies 
is defined as: 


40l | usually given in terms of the quantity S'(4) which 


S{J) = 


[E{L^) - E{L^ - 1 )] - [E{L^ - 1 ) - E{L^ - 2 )] 

mt) ' 


(5.33) 


where E stands for the absolute energy with respect to the ground state. In Ref. 


was 


shown that for triaxial 7 rigid cases 5'(4) > 0.56, with the limiting value corresponding to 
Z(4)-/?^ model. Studying Fig. |9l where ^(d) calculated in present model is visualized as 
function of a, one ascertains that the Z(4)-Sextic predictions fall in the aforementioned class 
for a > —2. Moreover, comparing present calculations with the value S{4) = 0.93 of the 
Z{4) solution one can see that it is doubly achievable in the a > —2 interval. In the rest of 
the a interval, ^(d) decreases to negative values as a —)■ — 00 , but not low enough to reach 
the U(5)-0(6) transition region values. The highest value of S'(4) obviously corresponds to 
ttc and which is very close to that of the Davydov’s triaxial rigid rotor model dlj. The 
phenomenon described above is known as the AJ = 1 or even-odd staggering. Taking 
another look at the a dependent spectra of Fig. |2l one can observe in the ground and 7 
bands another interesting phenomenon known as AJ = 2 staggering or AJ = 4 bifurcation 
which although very small was reported in the ground bands of actinide and rare earth nuclei 


42 


43| . There are many theoretical approaches dedicated to this topic which are briefly 


mentioned in Ref. 


43| . In the present model, this anomalous behaviour has a clear analytical 


origin which resides in the AL = 4 grouping of the states dehned by the rules (I3.16p and 
(I3.18p . It is interesting that the reciprocal closeness of the consecutive states is rearranged 
when going from negative to positive values of a. This theoretical result hints to the fact 
that the AJ = 2 staggering in the ground band of some nuclei can be due to higher order 
anharmonicities in their collective motion. 
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7.949 





FIG. 5. Energy spectra and some E2 transitions, normalized to the energy of the state 2^ and 
respectively to the transition probability B{E2,2'^ —0^), are visualized for each K = 1,2,3 and 
4 when a = Q. 


It is worth to mention that the similarities with the Z(4) and Z(4)-/3^ models enumerated 
so far reveal the fact that the approximation used to solve the eigenvalue problem for 
the 7 band is a good one. The advantage of the Z(4)-Sextic, comparing with the Z(4) 
and Z(4)-/3^ models, is that its potential can be varied smoothly, accommodating different 
deformation situations and creating in this way the possibility to cover intermediate cases 
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FIG. 6 . The same as in Fig. [H but for a = ac- 
between the Z(4)-/3^ and Z{4) or even beyond their bonndaries. 

TABLE I. Some B{E2) transitions, give n by Eq. (j4.81j) and normalized to B{E2;2^ —>• 0+), are 


compared with the experimental data 
predictions. 


45-47] for the I28,i30,i32^g isotopes and with the Z(4) model 


BiE2,Lt^L+) 

i28Xe 

i30Xe 

Z(4)-S 

i32xe 

Z(4) 

B{E2-,2+^0+) 

Exp. 

Z(4)-S^ 

Exp. 

Z(4)-S 

K ^ y 

1 

1 

1 

1 

1 

1 

^ 2 + 

1.468 

1.806 

1.966 

1.238 

2.048 

1.707 

69 + ^ 4+ 

1.940 

2.549 

2.972 


3.273 

2.414 

2 + ^ 2 + 

1.194 

1.771 

1.888 

1.775 

1.947 

1.737 

2 + ^ 0 + 

0.016 

0.000 

0.000 

0.003 

0.000 

0.000 


S is an abbreviation for Sextic. 


As in Ref. 


2l| . Z(4)-Sextic is applied in Fig. [7]for 128 , 130 , 132 ^^ isotopes which were consid¬ 


ered as candidates for Z(4) model and additionally in Fig. | 8 ]for I92,i94,i96p^ isotopes. The 
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FIG. 7. 
data 


'he theoretical energy spectra, given by Eq. (j3.27p . are compared with the experimental 


451447 1 of the I28,i30,i32^g isotopes and with the Z(4) model predictions. The corresponding 


rms values for Z(4)-Sextic are 0.516, 0.478 and 0.403, while for Z(4) one obtains 0.528, 0.389 and 
0.611. 


17 










































































































































12+ 9.769 

10 + 

8.624 

12 + 

9.253 

10+ 7.350 



10 + 

7.188 

8 + 5.773 

8 + 

6.377 

8 + 

5.324 

6 + 3.787 

6 + 

4.314 

6 + 

3.669 

4+ 2.439 

4+ 

2.479 

4"^ 

2.226 

2 + 1 

2 + 

1 

2 + 

1 

o 

+ 

o 

0 + 

0 

0 + 

0 

Z(4)-Sextic 
O' = 3.1 

Exp. 

Z(4) 


g 

band 




12 + 





10.116 

12 + 

9.253 


10 + 

8.672 


10 + 

7.573 





10 + 

7.188 


8 + 

6.392 


8 + 

5.880 





8 + 

5.324 


6 + 

4.298 


6 + 

3.835 

6 + 

3.669 

4+ 

2.415 4+ 

2.470.1 + 

2.226 

2 + 

1 2 + 

1 2 + 

1 

0 + 

0 0 + 

0 0 + 

0 

Z(4)-Sextic 

Exp. 

Z(4) 

a 

= 2.7 




g 

band 



12 + 

9.057 

10 + 

12 + 

8.558 

9.253 

10 + 

6.914 

8 + 

10 + 

6.333 

7.188 

8 + 

5.579 

8 + 

5.324 

6 + 

6 + 

3.709 

4.290 

6 + 

3.669 

4"^ 

2.513 4+ 

2.465 .,+ 

2.226 

2 + 

1 2 + 

1 2 + 

1 

0 + 

0 0 + 

0 0 + 

0 


Z(4)-Sextic Exp. Z(4) 
a = 43 

g band 


192pt 


8 + 

9.163 

8 + 

8.788 

9+ 

8.498 

9+ 

8.316 

6 + 

7^ 

6.358 — 

6.677 

QO'; 

6.357 

7+ 

6.065 

7+ 

6.188 

5+ 

5 + 

4.342 — 

4.682 

4.239 

4+ 

4.229 

3+ 

3.795 — 
4"^ 

2.910 , 

4.051 

3+ 

2.302 

3+ 

2.443 

2 + 

1.653 

1.935 1+ 

1.766 


Z(4)-Sextic Exp. Z(4) 
q' = 3.1 


7 band 


194pt 


7+ 

6.235 

7+ 

6.188 

5+ 

4.402 5+ 

4.563 5 , 

4.239 

4+ 

4.268 4+ 

3.743 — 
4+ 

4.051 


3+ 

2.809 ^ 


3+ 

2.332 


2.443 

2 + 

1.661 

1.894 0 + 

1.766 


Z(4)-Sextic Exp Z(4) 
a = 1.1 


yband 


196pt 


8 + 8.564 


9+ 7.943 
6 * 6.041 

7+ 5.737 
5 + 4.243 
4 + 4.179 

3 + 2.249 
? + 1.646 


8+ 

7.730 9 + 

8.316 


6 + 

6.357 

6 + 

5.644 7 + 

6.188 

5+ 

4.5265+ 

4.239 

4+ 

3.636^+ 

4.051 

3 + 

2.852 ^ 



3+ 

2.443 

2 + 

1.936 o+ 

1.766 


Z(4)-Sextic Exp Z(4) 
O'= 4.3 


yband 


4* 7.002 4^ 6.893 


2+ 4.995 
0+ 3.397 


2* 4.547 
0+ 3.776 


2+ 4.804 
0+ 2.954 


Z(4)-Sextic Exp. Z(4) 
Q' = 3.1 


ySband 


4+ 7.436 

4+ 6.893 


2+ 5.409 



2 + 

0+ 3.706 0+ 

4.603 

3.858 

2+ 4.804 

0+ 2.954 

Z(4)-Sextic 
a = 1.1 

Exp. 

Z(4) 


ySband 


4+ 6.893 

4+ 6.238 





2+ 4.804 

2+ 

4.308 




2+ 

3.828 


0+ 

2.954 0* 

3.192 

0+ 2.954 

z(4; 

I-Sextic 

Exp. 

Z(4) 

a 

= 4.3 




/3band 


FIG. 8. The same as in Fig. [71 but for the experimental data 


48-Mj of the I92,i94,i96p^ isotopes. 


The corresponding rms values for Z(4)-Sextic are 0.614, 0.543 and 0.682, while for Z(4) one obtains 
0.662, 0.707 and 0.732. 
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best fits of the experimental energy spectra were obtained for K = A, but not very different 
from K = 2 and 3, while the K = 1 case obviously has only theoretical importance. The 
experimental data of these nuclei are better described by Z(4)-Sextic comparing to Z(4), 
except for the ^^°Xe nucleus as can be deduced from the corresponding rms values given in 
the captions of Figs. [7] and [SI The ground and 7 bands of the Pt isotopes are similarly 
described by both models, while the (3 band is much better reproduced by Z(4)-Sextic. The 
poor description of the (3 band within Z (4) model is the reason why these nuclei were never 
considered as Z(4) candidates. In what concerns the Xe isotopes, their description by means 
of the present approach proved to be better than that of Z(4) only for ^^®Xe and ^^^Xe, 
although for ^^°Xe both models give very similar results. It is interesting that while the 7 
band of ^^®Xe has a similar description in both models, the ground and (3 bands are slightly 
better described within present approach. The picture is quite different in case of ^^^Xe 
where one obtains the biggest difference between the rms values, in favor of Z(4)-Sextic £t. 
Indeed, even if Z(4) describes the (3 band better it fails to do the same for ground and 7 
bands. Moreover, our model simulates very closely the AJ = 2 staggering of the experi- 



-10 -5 0 5 10 

a 

FIG. 9. The staggering 5'(4) given by Z(4)-Sextic for a G [—10,10] is compared with the values 
yielded by the Z(4) and Z(4)-/3^ models. 
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TABLE II. The same as in TableUbut for the experimental data 


48-50| of the I92,i94,i96p|. isotopes. 


B{E2,Lt ^L+) 

B{E2-,2+^0+) 


192pt 

194 

Pt 

196 

Pt 

Z(4) 

Exp. 

Z(4)-Sf. 

Exp. 

Z(4)-S 

Exp. 

Z(4)-S 

n ^ 0+ 

1 

1 

1 

1 

1 

1 

1 

^5 ^ 2 + 

1.556 

1.750 

1.728 

1.690 

1.478 

1.895 

1.706 

^9 ^ 4 + 

1.224 

2.424 

1.362 

2.296 

1.798 

2.770 

2.414 

8+ ^ 6 + 


3.078 

1.016 

2.880 

1.921 

3.622 

2.913 

lot ^ y 


3.493 

0.691 

3.238 


4.219 

3.293 

2+^0+ 


0.868 


0.810 

0.123 

1.047 

0.769 

2 ? ^ 7 


0.351 


0.275 

0.003 

0.590 

0.422 

^ OJ 


0.001 


0.001 

0.000 

0.000 

0.005 

2+ ^ 2t 


0.099 


0.088 

0.006 

0.126 

0.184 

Oj ^ 2 + 


1.362 

0.013 

1.159 

0.069 

1.852 

1.151 

0+^2+ 


0.000 

0.171 

0.000 

0.443 

0.000 

0.000 

6+ ^ 4+ 


1.028 


0.974 

1.207 

1.175 

1.142 

4+^2+ 


0.750 

0.427 

0.724 

0.714 

0.812 

0.801 

3+^2+ 

1.783 

2.251 


2.183 


2.415 

2.365 

6+ ^ 6 + 


0.224 


0.211 

0.394 

0.262 

0.218 

6+ — )• 4+ 

7 a 


0.000 


0.000 

0.012 

0.000 

0.000 

4+ 4+ 

7 ^ a 


0.365 

0.285 

0.347 


0.415 

0.381 

7 ^ 2+ 


0.000 

0.007 

0.000 

0.014 

0.000 

0.000 

3+ ^ 4+ 

0.664 

1.339 


1.280 


1.489 

1.360 

3+ ^ 2+ 

0.012 

0.000 


0.000 


0.000 

0.000 

2t ^ 2 * 


1.730 

1.809 

1.684 


1.837 

1.737 

2+ ^ 0 + 

0.010 

0.000 

0.006 

0.000 


0.000 

0.000 


S is an abbreviation for Sextic. 


mental ground band. The overall impression of the comparison with experiment is that the 
ground and (3 bands are quite well described for all considered nuclei, while the specihc 7 
band staggering of the rigid triaxial models in general, is found only in the experimental 
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spectrum of Plotting in Fig. [10] the potential fl3.20p using the htted values of a, one 

observes the following features. The advantage of Z(4)-Sextic over the pure oscillator or 
inhnite square well potentials, is that it can describe a shape phase transition. This can be 
clearly seen from the shapes of the potentials resulted for Xe isotopes, where ^^®Xe has a 
potential with a deformed minimum, ^^^Xe with a spherical minimum while ^^°Xe is situated 
in the critical point region. This result conhrms ^^°Xe as being a good candidate for the 
Z(4) model. Another important point emerging from Fig. [10] is that using a sextic potential 
with a deformed minimum one obtained a relatively good agreement with experiment for Pt 
isotopes, which could not be achieved using simpler potentials such as harmonic oscillator 
and inhnite square well. 



FIG. 10. Ground state potential (|3.20p with K = A plotted as function of y for all considered nuclei 
with a resulting from the fits of Figs. [7] and [8] The critical point potential (a = ac) as well as a 
potential close to the spherical limit (a = 20) are given for reference. 

Concerning the B{E2) transitions, the agreement of Z(4)-Sextic is very good with all the 
available data for the 128,130,132^^ isotopes. In what concerns the I92,i94,i96p|. isotopes, in 
the ground band the Z(4)-Sextic and Z(4) numerical results provide a good agrement with 
experiment only for the 4+ —>■ 2+ transition, the rest of the transition probabilities being 
overestimated in both calculations. A possible way to improve the agreement is to add 
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anharmonicities to the transition operator as was made in Ref. For transitions in the 
7 band and from the 7 band to the gronnd band both approaches give good resnlts, while 
for transitions from the (3 band to the gronnd and 7 bands the agreement is only partially 
good. These applications show that these isotopes can be considered partial candidates for 
Z(4)-Sextic and Z(4). The good agreement for all three bands of the isotope ^®^Pt proves 
that these solntions can describe real sitnations and opens the qnestion if there are better 
candidates. 


VI. CONCLUSIONS 


The main resnlt of the present work consists in the proposal of a new solntion for the 
Davydov-Chaban Hamiltonian, with a sextic oscillator potential for the variable [3 and 7 
’’frozen” to 30°. The solntion is conventionally called Z(4)-Sextic, in connection with the 
precedent Z(4) solntion where an inhnite sqnare well potential was considered. Choosing 
a qnasi-exactly solvable form for the sextic potential, a hnite set of states was analytically 
determined. The corresponding eigenvalne problem is exactly solved in the case of gronnd 
and (3 bands, while for the 7 band states an approximation is involved. The difference from 
the former qnasi-exactly solvable sextic potential approaches 28l-l31|. is the introdnction of 
a completely different scheme for angnlar momentnm attribntion which satisfy the condition 
of constant potential. Also it is the hrst time when the scaling property of the problem is 
employed to describe the properties of the qnasi-exactly solvable sextic potential. Indeed, 
as was shown in Section III, the model depends np to a scaling factor on a single parameter. 
Taking advantage of this property, one stndied the evolntion of the energy spectra and the 
corresponding transition probabilities when the free parameter is varied throngh different 
shapes of the associated sextic potential. For two valnes of the free parameter, the potential 
has one vanishing term. The spectra normalized to the energy of the hrst exited state and 
the B{E2) transitions normalized to the transition between the hrst excited state and the 
gronnd state calcnlated with the present model for these special cases, constitnte parameter 
independent realizations of the associated simplihed sextic potentials. 

A detailed comparison to the Z(4) and Z(4)-/5^ models, especially in terms of the energy 
spectrnm, revealed that the present formalism approximate qnite well the former in its 
critical point, and exactly reprodnces the latter in the asymptotic limit of the free parameter. 
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These facts suggest the consistency of the approximation used for the treatment of the 7 
band states. 

Numerical applications were performed for 128,130,132^^ I92,i94,i96p|- jgotopes. The 

results for the Xe isotopes revealed a shape phase transition with its critical point identihed 
with the ^^°Xe nucleus. The hts have a qualitative character, showing that the experimental 
realization of triaxial 7 rigidity is very much possible. Especially encouraging in this sense 
is the reproduction of the ^^^Pt spectrum. 

Concluding, one should say that the theoretical value of the proposed model resides in 
the fact that it adds to the few exactly solvable solutions of the collective model concerning 
only the ground and (3 bands, while its special cases contribute to the even more restrained 
set of parameter free models. 
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